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We present a non-perturbative method, called Parametric Perturbation Theory (PPT), which is 
alternative to the ordinary perturbation theory. The method relies on a principle of simplicity for 
the observable solutions, which are constrained to be linear in a certain (unphysical) parameter. 
The perturbative expansion is carried out in this parameter and not in the physical coupling (as 
in ordinary perturbation theory). We show that the method is capable to resum the divergent 
perturbative series, to extract the leading asymptotic (strong coupling) behavior and predict with 
high accuracy the coefficients of the perturbative series. 



In this letter we present a method, called Paramet- 
ric Perturbation Theory (PPT), which can be used to 
resum series, either divergent or with a finite radius of 
convergence, which appear in the perturbative solution 
of physical problems. The approach behind our method 
is completely new and it is based on few simple ideas: the 
first idea, which we call Principle of Absolute Simplicity 
(PAS), is that, instead of expanding the observable (en- 
ergy, frequency, . . .) in the physical coupling g, and thus 
obtain to a finite order a polynomial in g, we impose that 
the observable has the simplest possible form (linear) in 
a given unphysical parameter g; the second idea is that 
the expansion must be carried out in g and that the func- 
tional relation g = g(g) (unknown) must comply with the 
PAS to the order to which the calculation is done. This 
will allow to determine the relation between g and p and 
in turn to obtain the observable as a parametric function 
of g. 

Consider for example a model which depends on a pa- 
rameter g, and which is solvable when g = 0. Clearly, the 
application of Perturbation Theory (PT) to the calcula- 
tion of a physical observable O to a finite order yields a 
polynomial in g. Calling r the radius of convergence of 
the perturbative series, the direct use of PT must be re- 
stricted to \g\ < r. However, the misbehavior of the per- 
turbative series for O is the result of having expanded 
in a parameter, <?, which is not optimal. As a matter 
of fact, if one knew the exact solution to the problem, 
i.e. O = /(<?), then this solution could be considered 
as a polynomial of order one in the variable g = /(<?). 
Although this observation by itself cannot be used as a 
constructive principle, we may adopt the philosophy that 
the perturbative series for the observable can be simpler 
and convergent in all the domain, if it is cast in terms 
of a suitable parameter g. Only if such parameter, by 
luck or ability, turns out to be the g = f(g) discussed 
above, the exact solution is obtained. The goal, there- 
fore, is to progressively build this parameter g to yield 
an expression for O as simple as possible. In this frame- 
work the perturbative expansion is carried out in g and 
all the physical quantities in the problem are expressed as 
functions of g. In particular we have now that g = g(g). 



While the ordinary perturbation theory works by calcu- 
lating the contributions to higher orders in <?, each term 
of higher order refining the result to lower order, the ap- 
proach approach is the opposite: we carry out a per- 
turbative calculation in g, and then determine order by 
order the form of g = g(g) so that the observable O(g) 
can be a order one polynomial in g, as required by the 
Principle of Absolute Simplicity. 

Although in this letter we focus on the implementation 
of this method as a technique for resumming perturba- 
tive series, we show in a companion and more detailed 
paper that the same philosophy can be used to obtain a 
perturbation scheme fully autonomous from perturbation 
theory. 

Let us first sketch briefly how the method works and 
then apply it to some non-trivial perturbative series. 
Consider an observable O, represented through the per- 
turbative series perturbatively to some order 



A' 



(1) 



n=0 



where N = oo in some cases. The series could be either 
convergent, with a finite radius of convergence, or diver- 
gent, although we will not worry at this time. The full 
implementation of the method requires the introduction 
of an unphysical parameter, g, and the specification of a 
functional relation between g and g: 



9 = 9{q) 



(2) 



The choice of this relation is not completely arbitrary, 
since it must be capable of reproducing the perturbative 
terms when expanded around g = 0. For example, in 
many of the cases which we have studied we have used 



9{q) = Q 



1 + E 



N„ 



1 
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where the coefficients c n and d n are unknown to be de- 
termined applying the PAS. N u and Nd are integers. 

The Principle of Absolute Simplicity requires that the 
observable be linear in g, i.e. 



O^=bo + b ig . 



(4) 
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If wc substitute cq. ([3]) inside cq. JTJ) and expand 
around g = 0, working to a given order specified by the 
sum N u + Nd = N, then we can fully determine the co- 
efficients c n and d n by requiring that an equal number 
of nonlinear terms in g vanish, starting with the term 
going as g 2 . Notice that the choice of the integer param- 
eters N u and Nd determines the asymptotic behavior of 
g as g — ► oo (we are assuming for simplicity that the 
denominator in ([3]) never vanish for g > 0): 



g(o) ~ q 



N u -N d + l 



and in turn 



0< PPT ) « b + 6^1/(^.-^+1) 



(5) 



(0) 



Therefore in cases where the asymptotic behavior is 
known, one can choose N u and Nd to reproduce the exact 
asymptotic behavior of the solution; on the other hand, 
when the asymptotic behavior is unknown, working to 
a given order N u + Nd = N, one can select the most 
appropiatc asympotic behavior among those allowed by 
the combination of N u and Nd which keep N fixed. 

We will first apply the method to a model of (/> 4 field 
theory in zero dimensions , which will allow us to dis- 
cuss some other properties of our resummation. We con- 
sider the integral 



r+oo 

E(g) = / ■ 
Jo 



dx , 



(7) 



whose perturbative series is divergent. Eq.(|7| admits an 
exact analytical solution given by 



e 8 s / 1 

E( 9 ) = Wg K 1/4 



(8) 



where K\u(g) is the Bessel function of order 1/4. No- 
tice that for negative values of g this expression acquires 
an imaginary part, signaling that the system becomes 
metastable. 

We will now analyze this problem with the help of 
PPT, using a slightly different functional relation than 
the one in eq. ([3]): 



9{q) = Q 



1 



•£n=l Cne" 
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Notice that this relation implements the correct asymp- 
totic behavior, E(g) cx <? -1 / 4 as g — > oo. 

Using N — 3 we have found the transformation: 



1(Q) 



p (2925£ 3 + 881.8£> 2 + 58.8g + l) 5 
{-l7Z7.2g A + 2243.9g 3 + 832. 2g 2 + 58g + if 



Since go = —0.0593 is a zero of the denominator, 
lim e _,g g(g) = oo: this result signals the presence of 




FIG. 1: (color online) Comparison between the set [3,4] and 
the exact integral. The approximate solution has a branch 
point close to g — 0. 



a branch point in the proximity of go (see FigO}. The 
existence of this branch point can be understood as a sig- 
nal of the no-analyticity of the exact solution in g = 0, 
which is the reason why the perturbative series diverges. 

In the region g < the analytic continuation of the 
solution acquires an imaginary part. Using eq. fTT)]) we 
find the numerical solutions of g(g) = g, with g < 0. For 
example, corresponding to g = — 1 we obtain two pairs of 
complex conjugated roots accompanied by a single real 
root: 



qi = 0.20784 + 0.54897 i 

g 2 = 0.20784- 0.54897 ? 

Q 3 = -0.24838 

g 4 = -0.22321 + 0.00636i 

g 5 = -0.22321 - 0.00636i 



(11) 
(12) 

(13) 
(14) 
(15) 



Obtaining a complex value for g has an immediate ef- 
fect on the observable E(g), which acquires an imaginary 
part, ImE(g) = b\ Img. Therefore we can verify if one 
of these solutions corresponds to the exact solution of ([7]) 
for g = — 1: 



ImE{g) = -0.37679 



(16) 



which should be compared to the imaginary parts cal- 
culated with the PPT using the numerical roots gj, 
i= 1,...,5: 



EnE(g)r T 
ImE(g)P pT 
ImE{g)P pT 
ImE(g)r T 
ImE(g)r T 



-0.36488 
0.364886 


-0.00422 
0.00422 . 



(17) 
(18) 
(19) 
(20) 
(21) 
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TABLE I: Comparison between the perturbative coefficients 
of 0] and those predicted with PPT working with the set 
[3,2]. 





b 7 


bs 6 9 


bio 


611 


i exact 
On 


0.0fl06f39 


-0.00874935 0.00709675 


-0.00587143 


0.00493622 


6 [3,2] 


0.01106113 


-0.00874838 0.00709460 


-0.00586767 


0.00493037 


error (%) 


0.00233 


0.0ff08f739435 0.03022 


0.06403 


0.11851 




FIG. 2: (color online) Real and imaginary parts of E(g) for 
the zero dimensional 4 theory obtained using the PPT with 
N = 3. The results are compared with the exact expression 
of eq.(0. 



The comparison between exact and approximate real 
and imaginary parts of the integral for negative g is 
shown in FigO 

These results suggest that the first root corresponds 
to the analytic continuation of the solution for g > 
to negative values. On the other hand the real part of 
ReE{g) PPT = 0.74808 has the opposite sign ofReE(g) = 
—0.76033: this happens because our function is continous 
and therefore it is not possible to reproduce a disconti- 



nuity at g = 0. The exponential behavior of the exact 
solution for g — > 0~ cannot be reproduced in this ap- 
proach. 

We will now discuss a different problem, a lattice <fi 4 
model in 1 + 1 dimensions, described by the hamilto- 
nian Q 



2** 



15 



<t>i - 4>i+i ) + 7:' 



(22) 



where i is the site index and 7Tj and <f>i are canonically 
conjugated operators. 

Nishiyama has studied this model using both a linked 
cluster expansion, calculating the perturbative contribu- 
tions to order g 11 , and the Density Matrix Rcnormaliza- 
tion Group (DMRG). Since the perturbative series has 
a radius of convergence g f=a 1, he used an Aitken S 2 
process to accelerate the convergence of this series, ob- 
taining moderately improved results. Using these results, 
he speculated the existence of a singularity correspond- 
ing to g w —2 (using our notation) and of a Ising-type 
phase transition for a critical negative g. 

We can apply PPT to this problem using eq.© and 
considering the sets corresponding to N = 5. Comparing 
the difference E(g) PT — E(g) PPT we have found that the 
optimal set corresponds to \N U: N<j\ = [3,2]. In Table 
U we compare the exact perturbative coefficients going 
from 67 to 611 as given in [3j with those predicted by the 
set [3,2]. The largest error corresponds to bu and is of 
just 0.1%! 



In Fig|3]we have compared the perturbative polynomi- 
als for the energy from orders g 2 to g 11 with the energy 
resummed with the sets [3,2], [4,3] and [5,4]. There are 
several striking aspects which should impress the reader: 
first of all, the difference between the three sets is ex- 
tremely thin, thus signaling that the convergence is ex- 
tremely strong; in second place, the resummed energy 
confirms the DMRG result displayed in Fig. 2 of [§]; fi- 
nally, the resummed energy is a multivalued function, 
with a branch point at g w —1.025, falling extremely 



close to the singularity speculated by Nishiyama 1 . Be- 
cause of the use of a parameter g, PPT can deal with 
multivalued functions in a way which is not possible 
in conventional perturbation theory. Finally, the thiny 
dashed line in the plot corresponds to the numerical re- 
sult obtained in [3] using the Variational Sine Collocation 



1 Clearly, the branch point of a function y = f(x) at a point x = xq 
manifests itself as a singularity in that point when it is calculated 
using the Taylor series around a different point. 
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FIG. 3: (color online) Comparison between the resummed 
energies to orders [3,2], [4,3] and [5,4] and the perturbative 
polynomials. 



TABLE II: Virial coefficients for a hard spheres in 2 and 3 
dimensions given in Table I of 0] and predictions using PPT 
with different sets. 





B9/B2 

D 


B10/B2 
= 2 


Bg/ B2 

D 


B10/B2 
= 3 


Ref.[5j 


0.0362193 


0.0199537 


0.0013094 


0.0004035 


[0,5] 
[1,4] 
[2,3] 
[3,2] 
[4,1] 
[5,0] 


0.03739998 
0.03625994 
0.0362321 
0.0362551 
0.0368599 
0.1747048 


0.02496595 
0.02008503 
0.0199843 
0.02006717 
0.02258546 
0.75984885 


0.0023400 
0.0013509 
0.0013165 
0.0013404 
0.0017325 
0.0222648 


0.0031580 
0.0004884 
0.0004198 
0.0004664 
0.0014442 
0.0735264 



Method (VSCM) within a mean field approach. 

As a last example, we apply PPT to the prediction 
of virial coefficients of a hard sphere gas in D = 2, 3 
dimensions. Ref. [B| contains the values for the first 10 
virial coefficients for D = 2, . . . , 8. In this case, we 
have found out that the optimal choice corresponds to 
N u = Nd — 1, thus implying (using our notation) that 
linig^oo g(g) = g < 00, i.e. that the resummed function 



will have a singularity precisely at g = g. 
Working with the set [3, 4] we have found 

« 1.1625 , g%>%~ 1.43439. (23) 

In Table [II] we have considered the different sets cor- 
responding to N u + = 5, which use the first eigth 
virial coefficients as input, and used them to predict Bg 
and £> 10 , which have already been calculated [f|. The 
error (in percent) of the prediction of Bg with the set 
[2,3] is of about 0.035% and 0.54% in D = 2 and D = 3 
respectively. 

In Table IIIII we have compared our predictions for the 
virial coefficients B\\ through Bis using PPT with the 
set [3,4] with the predictions made in [5j, finding very 
similar results for D = 2 but rather different results in 
the case D = 3. 

Concluding, our method has several interesting fea- 
tures: it can describe multivalued functions, it can pro- 
vide the imaginary part of the observable corresponding 
to a mctastable state, it can resum divergent series and 
select the most appropriate asymptotic behavior of the 
solution among those available at a given order. Finally, 
it can also be used to make accurate predictions of yet 
unknown perturbative coefficients. Further applications 
of this method have been considered in a lengthier and 
more detailed paper. 
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TABLE III: Predicted coefficients for approximants with fO 
exact coefficients for D — 2 and D = 3. Comparison between 
the predictions of 0] and the predictions obtained using the 
set [3,4]. 





Bn/Bl° 


B12/S2 1 




B14/B2 3 


Bis/B 2 14 


Bie/B^ 5 


B1-/BI 6 


Sis/B 47 


D = 2 
Ref.[5] 
[3,4] 


1.089 x 10~ 2 
1.0901 x 10~ 2 


5.90 x 10~ 3 
5.9235 x 10~ 3 


3.18 x 10~ 3 
3.2117 x 10" 3 


1.70 x 10~ 3 
1.7421 x 10" 3 


9.10 x 10~ 4 
9.4698 x 10~ 4 


4.84 x 10~ 4 
5.1638 x 10~ 4 


2.56 x 10~ 4 
2.8247 x 10" 4 


1.36 x 10" 4 
1.5492 x 10~ 4 


D = 3 
Ref.[5] 
[3,4] 


1.22 x 10~ 4 
1.1599 x 10~ 4 


3.64 x 10~ 5 
2.2229 x 10~ 5 


1.08 x 10~ 5 
-8.5616 x 10~ 6 


3.2 x 10~ 6 
-1.8088 x 10~ 5 


9.2 x 10~ 7 
-2.0325 x 10~ 5 


2.6 x 10~ 7 
-2.0112 x 10~ 5 


-1.9136 x 10~ 5 


-1.7971 x 10~ 5 



